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Lyapunov-based quantum synchronization in a designed optomechanical system 
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We extend the concepts of quantum complete synchronization and phase synchronization, which 
are proposed firstly in [Phys. Rev. Lett, 111 103605 (2013)], to more widespread quantum gener¬ 
alized synchronization. The generalized synchronization can be considered as a necessary condition 
or a more flexible derivative of complete synchronization, and its criterion and synchronization 
measurement are further proposed and analyzed in this paper. As an example, we consider two 
typical generalized synchronizations in a designed optomechanical system. Unlike the effort to con¬ 
struct a special coupling synchronization system, we purposefully design extra control fields based 
on Lyapunov control theory. We find that the Lyapunov function can adapt to more flexible control 
objectives, which is more suitable for generalized synchronization control, and the control fields can 
be achieved simply with a time-variant voltage. Finally, the existence of quantum entanglement in 
different generalized synchronizations is also discussed. 

PACS numbers: 42.50.Wk, 05.45.Xt, 05.45.Mt, 03.65.Ud 


I. INTRODUCTION 

Complete synchronization and phase synchronization 
between two continuo us variable (CV) quantum systems 
were first studied W iMari et ^ in mesoscopic optome¬ 
chanical systems [j, and they also made the forward- 
looking prediction that the quantum synchronization is 
of potential but important applications in quantum in¬ 
formation processing (QIP). Subsequently, quantum syn¬ 
chronization is deeply discussed in variety of quantum 
systems, such as cav^ quantum electrodynamics [3, 
atomic ensembles [1, Q, van der Pol (VdP) oscillators 
, Bose-Einstein condensation Q ,superconducting 
circuit system Q, and so on. In these works, quantum 
synchronization is further extended from CV system to fi¬ 
nite dimensional Hilbert space corresponding to more ex¬ 
cellent quantum properties and the quantum correlation 
is also analyzed quantificationally in those synchronous 
quantum system In addition, quantum syn¬ 

chronization criteria U 0, HI, fl^ and synchronization 
between nodes in quantum network are still hot topics in 
the field of quantum synchronization theory. 

Generally, existing quantum synchronization schemes 
can be attributed to the concept of coupling synchro¬ 
nization, i.e., one subsystem of the synchronous systems 
plays the role of controller acting on the other subsys¬ 
tem [IH3,0, ini- A significant advantage of this kind 
of direct linking is its high maneuverability. However, it 
still remains some difficulties to achieve the better appli¬ 
cations in QIP with quantum synchronization. For the 
weak coupling of the quantum level, it is difficult to elim¬ 
inate the difference between systems if it is big enough, 
and in fact, it is often in this case. Fundamentally, too 
strong driving or pump fields will compel systems to the 
form of the forced synchronization, just like what they 
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are realized in previous works [MO- This deficiency 
of coupling synchronization is a severe limitation, which 
causes other types of synchronizations are hardly realized 
in addition to complete and phase synchronizations. For 
examples, antiphase synchronization and projective syn¬ 
chronization, which are also widely applied in the classi¬ 
cal synchronization fields [l3l - [T5l| . were rarely discussed 
in quantum systems. 

In traditional control theory, besides the coupling 
terms, there exists an external controller which is im¬ 
posed on response system in order to provide a more 
outstanding control capability, implying that a designed 
controller can establish a more flexible relationship be¬ 
tween two controlled subsystems [l^ . It enlightens us to 
think about such problems: can a more generalized syn¬ 
chronization (like above mentioned antiphase and pro¬ 
jective synchronizations) be extended and obtained in 
quantum domain? If it works, what kinds of the crite¬ 
rion and measurement are needed in this quantum gen¬ 
eralized synchronization? And most importantly, how 
are the controllers designed in order to satisfy various 
requirements corresponding to different kinds of general¬ 
ized synchronizations? 

For responsing above questions, in this paper, we study 
the general properties of different synchronization forms 
and expand them into quantum mechanics based on 
Mari’s complete synchronization theory. The criterion 
and measurement of the generalized synchronization are 
also proposed and they will be divided into two orders 
for calculating and analyzing conveniently. Instead of di¬ 
rectly establishing interaction between two subsystems, 
here we utilize Lyapunov control theory which has exhib¬ 
ited comprehensive applications in target quantum state 
preparation and suppressing decoherence to design the 
external controller [la-llq . Although the Lyapunov func¬ 
tion is constituted by expectation values, our results show 
that the quantum fluctuations can also be effectively sub¬ 
dued by the controller. In addition, the classical and 
quantum correlations are considered via calculating the 
Lyapunov exponent and Gaussian Negativity. In par- 
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ticular, we demonstrate that CV entanglement can exist 
in generalized synchronization, however, it will disappear 
if generalized synchronization tends to complete synchro¬ 
nization. This phenomenon is consistent with Mari’s and 
Ameri’s conclusions about entanglement in complete syn¬ 
chronization Qi- Therefore, we believe existing quan¬ 
tum complete synchronization can be included in our gen¬ 
eralized synchronization theory. 

We organize this paper as follows: in Sec. m we in¬ 
troduce the definition and the properties, especially mea¬ 
surement method, of quantum generalized synchroniza¬ 
tion. In Sec. mi we analyze the dynamics of an optome¬ 
chanical system, and realize two kinds of representative 
generalized synchronization (constant error synchroniza¬ 
tion in Sec. IIII Al and time delay synchronization in Sec. 
IIII B|) respectively via designing appropriate control field 
based on Lyapunov function theory. The correlation in 
generalized synchronization is also discussed in Sec. IIVI 
and a summary is finally given in Sec. |Vl 


II. QUANTUM GENERALIZED 
SYNCHRONIZATION 

We begin this section with a brief introduction of gen¬ 
eralized synchronization and its expansion in quantum 
domain. Considering two general classical dynamics sys¬ 
tems whose evolutions satisfy following equations: 

dtXi{t) = F{xi{t)) + Ucl{xi,X2) + Uel 
dtX2it) = F{x2(t)) + Uc2{xi,X2) + Ue2 

here 3:1,2 (t) G are the state variables of two systems in 
time t. Uci,2 are the mutual couplings between systems 
and correspondingly, C/ei,2 are the external controllers 
belonging to their respective system. If there are con¬ 
tinuous mappings /ii, /12 : R^ —>■ R^ and following syn¬ 
chronization condition in Eq. © can be achieved when 
t —>■ 00, then two systems depending on /ii, /12, and 
X2 will be of the consistent evolution. 

lim \hi{xi{t)) - h2{x2{t))\ 0 (2) 

This controllable correlation is named as generalized syn¬ 
chronization, and it will degenerate to common complete 
synchronization or phase synchronization via selecting 
hi(xi) = Xi or hi{xi) = axg{xi), respectively. Simila^ 
with Mari’s measurement Sc{t) := {qtit) + [l[, 

generalized synchronization can be extended from clas¬ 
sical to quantum by considering conjugate quantities si¬ 
multaneously, and the corresponding measurement can 
be defined as: 

( 3 ) 

where gg_ := (/ii(gi) -/i2(g2))/A/2 and Pg_ := (/ii(pi)- 
^2(P2))/\/2 are the quantized generalized error operators. 

Nevertheless, it is not easy to use Eq. ([ 3 |) directly in a 


concrete model. In some cases, /ii,2(9,Pi,2) are not strict 
physical descriptions because they are actually superop¬ 
erators. On the other hand, it could be difficult to calcu¬ 
late Sg{t) in CV quantum systems. Therefore, in order 
to analyze quantum synchronization in CV mesoscopic 
systems, we adopt mean-field approximation to simplify 
Eq. (j 31 ). Then synchronization measurement can be di¬ 
vided into two parts: the first order criteria is to describe 
the consistency of expectation values: 

lim \hi{qi{t)) - /i 2 ( 92 (t))| 0 
lim \hi{pi{t)) - h2ip2{t))\ 0 

t—^OO 

and the second order measurement is to determine fol¬ 
lowing quantum fluctuations: 

S'it):={6ql_{t) + Spl_{t))-^ ( 5 ) 

where o refers to (o) and So := 0 — o ior o € {qg-,Pg-}. 

Physical meaning of Eq. ([ 4 ]) and Eq. ([ 5 ]) are more def¬ 
inite to explain quantum synchronization, i.e., systems’ 
expectation values are required to satisfy the “classical” 
generalized synchronization conditions and the pertur¬ 
bation on synchronization behavior caused by quantum 
effect is squeezed as much as possible. To verify this, 
Eq. ® will be equivalent to Eq. ([ 3 ]) if the first order 
criteria are satisfied. Conversely, if a designed external 
field can not only make the evolution of systems to re¬ 
alize “classical” generalized synchronization conditions, 
but also increase corresponding second order measure¬ 
ment Sg, then it can be seen as an appropriate control 
field for realizing quantum synchronization. This is the 
basic idea of designing the control field. 

In some particular models, if hi and /12 are selected 
as flat mappings, Eq. can be further simplified to 
measure fluctuation 

S^git):={SqUt) + 5pl{t))-^ ( 6 ) 

Here g_ := {qi — q2)lV^ and p_ := (pi —p2)lV^- Com¬ 
pared with Eq. ([ 5 ]), Eq. ([ 3 ]) is more easy to be obtained 
via the covariance matrix of system. 


III. LYAPUNOV-BASED SYNCHRONIZATION 
IN OPTOMECHANICAL SYSTEM 

We analyze Lyapunov-based synchronization in op¬ 
tomechanical system to more intuitively explain above 
theory of quantum generalized synchronization. Our 
model consists of two oscillators which couple with a 
Fabry-Perot cavity together (see Fig. [Ij. The Hamil¬ 
tonian corresponding to this model can be divide into 
four parts: F[ = Hq + Flmt + H^v + Hdt). Here 
Ho = ojia^a + is a sum of free 

Hamiltonians corresponding to the optical field and two 
oscillators. Moreover, Hint = —gia^aqi — 520^0(72 and 
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FIG. 1. Diagram of optomechanical system corresponding to 
our model. Here two oscillators are placed at wave nodes of 
a Fabry-Perot cavity and they couple with the cavity field 
via linear optomechanical interactions, and their origins are 
respectively set at the equilibrium positions. 


Hdiv = iE{a ^are the standard forms of 
optomechanical interaction and driving field respectively 
(l 9 l| . Hc{t) is an external control Hamiltonian which rep¬ 
resents the coupling with designed time-dependent fields. 
In our model, we consider such a form of control field 
which can create a deviation in respective potential term 
of two oscillators. This effect can be regarded as a time- 
dependent rescaling of the mirror frequency [l^, i.e., 

+ (7) 

We will give a more detailed discussion about how to 
realize this form of control field using specific experiments 
in Sec. El Then the Hamiltonian of the whole system is 
written as follow after a frame rotating: 

H=J2 + ^[1 + C,[t)]q^^ - gjO^aqj} 

i=i.2 ( 8 ) 

— Aa^a -b iE{a^ — a) 

In above expressions, a(a^) is the annihilation (creation) 
operator for the optical field and correspondingly for 
j = 1 , 2 , qj and pj are dimensionless position and momen¬ 
tum operators of the oscillator j respectively. A = w; 
refers to the detuning between the frequencies of the laser 
drive and the cavity mode, ojmj is the mechanical fre¬ 
quency. Pj is the optomechanical coupling constant and 
E is the drive intensity. In order to solve the dynam¬ 
ics of the system, we consider the dissipative effects in 
the Heisenberg picture and write the quantum Langevin 
equations as follows [ 20 l - l^ : 

dta = (—K + iA)a + igiaqi + *52092 + E + 

^tqj — ^mjPj ( 9 ) 

dtPj = + Cj(t)]qj - jjPj + gjtt^a + ij 

Here k is the decay rate of optical cavity, and 7^ is the 
mechanical damping rate of each oscillator, a™ is the 
input bath operator, which satisfies (a®"(t)a™’^(t')) = 
5 {t — t') [2^. Similarly, ^j{t) is the Brownian noise op¬ 


erator describing the dissipative friction force acting on 
the jth mirror. In the Markovian approximation, the 
autocorrelation function of ^j(t) satisfies the relation: 

= 7^(2h6 -b l)5jfS{t - t'), 
where Uf, = [exp^fiLUj/ksT) — 1 ]“^ is the mean phonon 
number of the mechanical bath which gauges the tem¬ 
perature T [234^ . 

Solving directly a set of nonlinear differential operator 
equations like Eq. Q is quite difficult, however, a mean- 
field approximation is a ccep table in our mesoscopic op¬ 
tomechanical model l 2 u - l^ . On the other hand, as 
we discussed in Sec. im the quantum synchronization 
measurement modifying by mean-field approximation 
can describe the generalized synchronization effect more 
accurately. Therefore, the every operator in Eq. m can 
be rewritten respectively as a sum of its expectation value 
and a small fluctuation near the expectation value, that 
is, a{t) = A{t) +Sa{t), d{t) = o{t) +So{t),o S (91.2,Pi,2)- 
After neglecting the high-order fluctuation terms, the 
“classical” properties of our optomechanical system can 
be described by following nonlinear equations: 

dtA = {—K + iA)A + igiAqi + ig2Aq2 + E 

^tPj — ^mjPj (fO) 

dtPj = + Cj{t)]qj - jiPj -b gj\A\‘^ 

and the corresponding quantum fluctuations can also be 
confirmed by: 

dtSa ={—K + iA)Sa + igj{qjSa + ASqj) + 

i=i,2 

dtSqj =uJmjSpj 

dtSpj = - + Cj{t)]Sqj - "ijSpj -b gj{A*6a + A5a)) -b ij 

( 11 ) 

Transforming the annihilation operators as the forms 
of a = (x -b iy)ly/2 and a™ = (i*" -b *5“)/V2 

respectively. Then Eq. m can be rewritten 
more concisely as dtu = Su -b C by setting the 
vectors u = {Sx,Sy,Sqi,6pi,6q2,Sp2)~^ and i = 

(x™, 9™, 0,^1, 0 , ^2)^, and the corresponding S' is a time- 
dependent coefficient matrix (see Appendix. A for more 
details). In this representation, the evolution of correla¬ 
tion matrix D defined as 

Dij{t) = Dj^{t) = ]^{ui{t)Uj{t) +Uj{t)Ur{t)) (12) 

can be derived directly by (see [illiJii) 

dtD = SD + DS^ -b N (13) 

N in Eq. (ITSl) is a noise matrix and it will be a diagonal 
form, i.e., diag(«:, k, 0,7i(2hb -b l),0,72(2hb -b 1)), if the 
noise correlation is defined by (Ci(^)Ci(l) + Cj(l)Ci(0)/2 = 
Nij5{t — t'). With the help of Eq. (|T^ . above mentioned 
synchronization measurement S' can be simply expressed 
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as 


S'^it) ={6q^_{t)+SpUt))-^ 

={2 [^ 33(0 + D55{t) - 2D35{t)] ('X4) 

+—[D 44 {t) + DQQ{t) — 2 D 4 Q{t)]} ^ 

and its evolution can be obtained by solving Eq. (nni) and 
Eq. (I13|) in order. At this point, all dynamic properties 
of our system, including synchronization and correlation, 
can be learned by means of the solutions of Eq. 
m and (fT^ . In the following subsections, we introduce 
two common forms of generalized synchronization: con¬ 
stant error synchronization (IIII Al) and time delay syn¬ 
chronization (IIII B|) to exhibit our ability about control 
synchronization, we will also prove how the controller is 
designed to realize these synchronizations. 


A. Constant error synchronization 


Constant error synchronization can be regarded as a 
translation in phase space between two systems. In Eq. 
m, if we let hi{xi) = xi + ci and 112(3:2) = X2 + C2, then 
the “classical” synchronization criterion will be 

lim |a;i(t) - a;2(f)|-)■ C2 - Cl = c_ (15) 

where c_ is the so-called constant error. In view of the 
controller influencing directly on dtPj , we firstly only con¬ 
sider the evolutions of the momentum operators, and fur¬ 
ther construct following Lyapunov function by using their 
expectation values: 

Vp{t) = {pi{t) - p 2 {t)f ( 16 ) 

One can easily verify that Vp 0 and V), = 0 is valid 
only when pi{t) — P 2 {t) = 0 . Substituting Eq. (ITUl) into 
Eg. lfTHl) . the time derivative of Vp can be calculated hand¬ 
ily if Ci{t) = C2{t) = C{t) 

Vp{t) =2{pi{t) - P2{t)){pi{t) - P2it)) 

=2{[1 + C{t)]{uJm2q2 - UJmiqi) “ 7lPl (17) 
+ I2P2 + (51 -52)|Ap}(pi -P2) 

we find that Vp{t) is always non-positive by setting 

pi{t)-p2{t) =-P2{t)) ( 18 ) 

where fc is a positive real number. With this choice, Vp 
simultaneously satisfies Vp ^ 0 and Vp = — 2 k{pi{t) — 
P 2 {t)Y ^ 0 - Under this condition, the system will grad¬ 
ually evolve to a stable state which corresponds to the 
origin of the Lyapunov function, i.e, pi{t) = P2{t) [ 3 - 
In order to satisfy the required form of Lyapunov func¬ 
tion, the control field can be obtained based on Eqs. (EZl) 


and (flSjl 


C{t) 


{l-k)[pi -P2] - (gi - g2)|Ap 

Wm 2(?2 — i^miqi 


(19) 


here we have already set 71 = 72 = 7- It is needed to 
emphasize that all mechanical quantities without specif¬ 
ically being marked in this equation represent the expec¬ 
tation values at time t (e.g., pi := Pi{t)). Otherwise, 
confusions may occur in following discussion about time 
delay synchronization. 

We notice, however, the control field C[t) in Eq. (fTOl) 
could be infinite when the tracks of qi and q2 are adja¬ 
cent. In particular, complete synchronization is not ac¬ 
ceptable if ojmi — ^m2- For avoiding this singularity, it 
is necessary to add a lower bound in the denominator of 
the control field. Therefore, the control field is modified 
in follow form: 


' (7 - k)[pi -P2] - {gi -52)|A|" 


- 1 


Wm 292 — ^miqi 

C{t) = {when \wm2q2 - Wmiqi \ > c_) ( 20 ) 

0 

^ {when \wm2q2 - Wmigil ^ c_) 


The physical mechanism corresponding to Eq. (EOl) can 
be interpreted as follows: Assuming the gap between two 
oscillators is small enough to satisfy \uJm2q2—Wmiqi \ ^ c_ 
at initial moment, then the control field will not work and 
the difference between oscillators under different Hamil- 



FIG. 2. Evolutions of expectation values (blue dashed and 
red solid), control field and errors (green dashed and black 
solid respectively denote control field is imposed or removed), 
in which (a) and (b) respectively corresponds to momentum 
and position operators of each oscillator. Here we set A = 1 
as a unit and other parameters are: uimi = 1, u)m2 = 1.005, 
g\ — 0.008, g2 = 0.005, E = IQ, k — 0.15, 7 = 0.005 and 
fib = 0.05. For the control field, the parameters are taken as 
k = 2 and C- = 3. (c): The limit cycles of two oscillators in 
phase space, (d): Robustness of our control system. 
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tonians will increase; Once such difference crosses the 
boundary \u!m2<l2 — t^mi 9 i| > c_, then non-zero control 
field will drag their orbits to close each other until a crit¬ 
ical distance is reached, which will cause the invalid con¬ 
trol field is resumed again. With time going by, the error 
evolution will be controlled in a stable limit ellipse. Un¬ 
der particular k and c_, it can be regard as a fixed point 
if the major axis of this ellipsoid is small enough. In this 
case, two systems will finally realize such a synchroniza¬ 
tion: Pi — P2 = 0 and — 52 = c_. Therefore, we make 
sure that C{t) in Eq. (ITOl) is able to control the system 
achieving generalized synchronization. 

In Fig. [ 5 ] we provide simulation results of the two os¬ 
cillators to verify the synchronization phenomenon un¬ 
der the control field. In Fig. (DJa), one can directly 
see momentums of two oscillators will take on consis¬ 
tent evolution after t = 41 . 3 , which is exactly the same 
with the time point that control field is non-zero. Corre¬ 
spondingly, the momentum error will stabilize at zero in¬ 
stead of generally enlarging along with the control field. 
Fig. [IJa) also quantitatively shows the control field is 
a slowly varying function of the time. Such a slowly 
varing control field can improve the stability of the sys¬ 
tem, simultaneously, it is more easily to be implemented 
by experiments. In Fig. [ 2 ljb), we plot the positions of 
two oscillators and corresponding errors. It illustrates 
that, although two oscillators are not consistent in their 
positions, the error can still maintain a constant (c_). 
Taken Fig. [SJa) and Fig. [ 2 jb) together, we can deter¬ 
mine constant error synchronization between two oscil¬ 
lators is achieved. In Fig. [H^c), we show the “tracks” 
of two oscillators in phase space. Two oscillators will 
evolve to their respective limit cycle and, as we pre¬ 
dicted above, constant error synchronization corresponds 
to a translation between limit cycles in phase space. Fig. 
[ 2 ]Jd) reports the robustness of our synchronization sys¬ 
tem. Here we assume each quantity in Eq. dsni) has 
been added a Gaussian noise whose standard deviation 
is a, i.e., o(t) = J\f{o(t),a) (o G {51^2,Pi,2,and when 
the final control field is added on the system, it also has 
a noise (C(t) = JV(C(t), a)). The accuracy of the syn¬ 
chronization scheme in this case is described by following 
auxiliary quantity 



FIG. 3. (a): Evolution of modified synchronization measure¬ 
ment. (b): Time-averaged synchronization measurement with 
varied bath temperature. Here blue of corresponds to control 
field imposed and red is the case control field disappears. 

creasing Sg substitutes the trend to 0 , which significantly 
outperforms uncontrolled situation. Therefore, we rec¬ 
ognize that the control field can indeed achieve quantum 
control rather than the synchronization of the classical 
level. In Fig. [ 31 [b), we also show how the bath temper¬ 
ature will influence synchronization phenomenon and it 
can be known that Sg(t) will keep almost unchanging if 
the bath temperature is limited within ImK^fib = 0.28 
corresponding to a MHz phonon frequency), and it is 
still larger than that belonging to the uncontrolled sys¬ 
tem even though T goes up to 10 mK{nb = 6 . 14 ). This 
range is qnite broad compared to other correlation con¬ 
trol schemes in optomechanical systems SEl. 

B. Time delay synchronization 

Time delay synchronization can be regarded as a con¬ 
stant phase deviation between two systems, and their 
“tracks” in phase space are overlapping like complete syn¬ 
chronization. In Eq. m, if we set hi{xi) = xi{t) and 
^2(3^2) = X2{t — t), then the “classical” synchronization 
criterion will be 

lim |a:i(f) - a: 2 (i - t)|-)> 0 ( 22 ) 

t—^OO ^ > 

Similarly with the above discussion, we define following 
Lyapunov function: 




where p'_,q'_ refer to the errors in biased control field, 
and r is the average radius of the limit cycle. One can 
find that R{a) will always remain above 96 % even the 
a — 0 . 02 . Under those parameters, even if there are ob¬ 
vious fluctuations in the control field, however, the errors 
between two oscillators are still stable to approach 0 , c'_. 
Therefore, we confirm our control is stable enough for 
some interferences. 

Besides the expectation values, Mari’s measurement is 
also calculated to prove that quantum fluctuation is simi¬ 
larly squeezed by control field. Fig. [ 3 Da) illustrates an in- 


Vpi.*) = - P 2 {t - T)f ( 23 ) 

and its derivative can also be expressed as I^ = 2 (pi — 
hit - T))iPi -P2{t - r)), where 

Pi-hit -t) = -Wmi[l + Ci]qi - 7 ipi -b gi\A\'^ 

+ Wm252(< -t)+ j2P2it - t) + 52|^(i “ t)^ 

It needs to emphasize again that, in above expres¬ 
sions, all mechanical quantities without specifically be¬ 
ing marked represent the expectation values at time t. 
Similarly, we let 

h-h{t -t) = -k{pi{t) - p 2 it - t)) ( 25 ) 
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to satisfy Vp ^ 0 and Vp = — 2 k(pi — p2(t — t))^ ^ 0 . 

I 


Then the corresponding control field will become: 


Ciit) 


(7 - k)[pi -P2{t - t )] - gi|Ap + g2\A{t - t)\^ - UJm2q2{t - t) 

—Wml?! 


( 26 ) 


via setting (72 (t) = 0 and 71 = 72 = 7 for simplicity. 

Eq. () 26 l) is also of singular point at qi(t) = 0 , therefore, an artificial boundary is necessary to avoid an infinite 
control field too. Unlike Eq. O, our purpose here is to make two systems achieve complete synchronization after 
eliminating the time delay. Therefore, this limitation is on the whole control field instead of the denominator. Then 
the control field should be 


Ci(t) 


(7 - fc)[Pi -P2(t - t )] - gi\A\'^ + g2\A{t - r)p - LUm2q2{t - t) 

—Wml9l 

Cm 

-Cm 


{—Cm < (7i ^ Cm) 


{Cl > Cm) 
{Cl < -Cm) 


( 27 ) 


In Fig. IlKa) and (b), we show that the evolution of one 
oscillator seems to be a time translation of the other oscil¬ 
lator, and the errors tend to zero like complete synchro¬ 
nization after eliminating the time delay. It also exhibits 
a quickly varying control field, which is different with 
the performance in constant error synchronization. Gen¬ 
erally speaking, quickly varying control field can make 
the system achieve synchronization faster. Fig. HKa) and 
(b) show two oscillators will achieve synchronization in 
a short period of time t < 10 , which shortens four times 
time relative to Fig. [ 5 ] (a) and (b). Furthermore, Fig. 
IDJd) illustrates synchronization accuracy is 96 %, which 
means robustness is also remained in a high level under 
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FIG. 4. Evolutions of expectation values (blue dashed and 
red solid), control field and errors (green dashed and black 
solid respectively denote control field is imposed or removed), 
in which (a) and (b) respectively corresponds to momentum 
and position operators of each oscillator, (c): The limit cycles 
of two oscillators in phase space, (d): Robustness of control 
system. Here we set r = 5, Cm = 1 and other parameters are 
the same with Fig. [3 


the fast oscillating control field. We also plot the limit 
cycles of two oscillators in Fig. SKc). It can be seen that 
two limit cycles are almost coincident in most of time 
except the short time intervals at origin and destination 
points in which both limit circles take on inconsistent 
evolutions because of the time delay. 

We also consider the evolution of quantum fluctuation. 
Fig. ISDa) shows that the fluctuations are further squeezed 
by quickly varying control field, and Sg{t) exhibits a 
higher value than that in constant error synchronization. 
Fig. [SKb) also illustrates that destruction on the synchro¬ 
nization effect causing by environment is also weakened 
and Sg{t) will still remain high even at T = lOmK. From 
this perspective, we believe that quickly varying control 
field is a appropriate form of synchronization control too. 


IV. CORRELATION IN GENERALIZED 
SYNCHRONIZATION 


The correlation between synchronized quantum sys¬ 
tems is a important research object in QIP. Intuitively, 
two different systems can achieve consistency in some ex- 




FIG. 5. (a): Evolution of modified synchronization measure¬ 
ment. (b): Time-averaged synchronization measurement with 
varied bath temperature. Here blue corresponds to control 
field imposed and red is the case control field disappears. 
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FIG. 6. The largest Lyapunov exponents of the errors with 
varied c_ (a), r and Cm (b) which corresponds to constant 
error synchronization and time delay synchronization, respec¬ 
tively. 




FIG. 7. The maximum Negativity of the errors with varied 
c_ (a), r and Cm (b) which corresponds to constant error 
synchronization and time delay synchronization, respectively. 


tent, meaning that there inevitably exists a certain cor¬ 
relation between those systems. To verify this, quantum 
mutual information which is a measurement of total cor¬ 
relation has been proved to be homology with synchro¬ 
nization measurement in VdP oscillators 0 . However, 
it is very difficult to identify the type of this correla¬ 
tion. Especially, whether quantum entanglement exists 
in the CV synchronization is controversial 0 , 11113 . 
Therefore, we pay more attention to the properties of en¬ 
tanglement when we consider the quantum correlation in 
our model. 

The mean-field approximation used above can make 
it more convenient to analyze classical correlation and 
quantum entanglement. The classical correlation can be 
verified via calculating the largest Lyapunov exponent of 
the errors [ll|, [ 13 , i.e., = max{Lj^(pg_), Ly ((jg_)}, 

where 


Lj,(o) = lim — In 


So{t) 

5 o( 0 ) 


{o & {Pg-,qg-}) ( 28 ) 


Correspondingly, CV quantum entanglement is measured 
by Gaussian Negativity = max{ 0 , — log2 [H, S, 
[33, here 
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A(r) = det A -I- det B — 2 det C and 
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systems. This means, compared to the direct coupling 
between two systems, designed Lyapunov function can 
more effectively help the system to establish correlation 
in the expect value level. 

In Fig. [3 we plot the results of maximum Negativity 
under different characteristic parameters. It should be 
noted that generalized synchronization is actually a nec¬ 
essary condition for complete synchronization and char¬ 
acteristic parameters can be regarded as a description 
about the gap between generalized synchronization and 
complete synchronization. Consequently, one can com¬ 
pare complete synchronization with generalized synchro¬ 
nization via setting c_ = 0 , r = 0 and Cm —t 00. Fig. 
Ela) shows that quantum entanglement does not exist 
until c_ = 2 . 65 . Then it takes on a rising trend and 
its the maximum value is at c_ = 3 . 13 . Subsequently, 
Negativity will decline and tend to zero again. The evo¬ 
lution of maximum Negativity can be understood as fol¬ 
low: when c_ is small enough, the oscillators are equiva¬ 
lent to achieving complete synchronization and they are 
always separable in this case. This opinion is the same 
with Mari’s conclusion, i.e., the evolution in superposi¬ 
tion tracks will hinder the generation of CV entangle¬ 
ment. With gradually going away of two tracks, entan¬ 
glement is allowed to exist because the oscillators are no 
longer complete synchronization. Finally, if the distance 
of two tracks increases constantly, quantum correlation 
between two oscillators is remarkably weakened and en¬ 
tanglement will disappeare again when c_ is large. In 
Fig. EDb), similar conclusions are also obtained in time 
delay synchronization. Therefore, we think that Gaus¬ 
sian entanglement can coexist with generalized synchro¬ 
nization, although it may be prohibited by complete syn¬ 
chronization. 


In Fig. El we plot the largest Lyapunov exponent under 
different characteristic parameters (c_, t and Cm)- The 
results show that negative Lyapunov exponent is always 
present under different control fields. This performance is 
superior to our conclusion in Ref. , in which we found 
that the positive-negative of Lyapunov exponent corre¬ 
sponding to coupling synchronization depends sensitively 
on varing coupling intensity between two opomechanical 


V. DISCUSSION AND RESULTS 

Here we give a brief discussion about the parameters 
of our optpmechanical system and the realization scheme 
of the control field. The parameters selected in the simu¬ 
lation are similar with Ref. 0, |l3,113, [13 113 • However, 
in order to highlight the roles of the coupling and the 






















































controller, we^propriately reduce the value of the driv¬ 
ing intensity [ll|. Beyond that, the deviation in poten¬ 
tial like Eg. o has been investigated by theoretical re¬ 
searches and recent works report that the deviation 
can be achieved by us ing charged m echanical resonators 
For example, IZhang et al.l found that the effec¬ 
tive frequency = ^^[(1 + <11(0)] can be controlled 
by a time-dependent bias gate voltage U{t) = Uof{t). 
And the relationship between dimensionless factor f{t) 
and C{t) is C{t) = T]f{t), where 


CqUqQmr 


(31) 


is obtained in Ref. [s^. Therefore, we are sure the con¬ 
trol terms in Eq. (EOl) and (l?71) cab be achieved easily in 
a specific experiment. 

In summary, we have extend Mari’s theories about 
quantum complete synchronization and phase synchro¬ 
nization to a more general situation that we defined 
as quantum generalized synchronization in this paper. 
The corresponding control methods, criteria and mea¬ 
surements are also proposed quantificationally based on 
Lyapunov function, Lyapunov exponent and modified 
Mari’s measurement. This generalized synchronization 
can be regard as a prerequisite of traditional quantum 
synchronization, and it can establish more flexible re¬ 
lations between two controlled systems. To verify this, 
we have shown that some important properties in our 


model, such as entanglement in synchronization, will be 
consistent with previous works if the generalized synchro¬ 
nization tends to complete synchronization. So, design¬ 
ers can complete different synchronizations according to 
their requirements based on our theory. For making our 
theory more intuitive, we have considered two common 
generalized synchronization, that is, so-called constant 
error synchronization and time delay synchronization in 
an optomechanical system. With the help of control fields 
designed by Lyapunov function, we have proved two os¬ 
cillators can satisfy the requirements of various synchro¬ 
nizations. We believe that our work can bring certain 
application values in quantum information transmission, 
quantum control, and quantum logical processing. 
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Appendix A: Parameters in quantum Langevin 
equations 

The concrete form of coefficient matrix S in Eq. m 
is 


S = 


f —K 
A + giqi + 3292 
0 

'j2giKe{A) 

0 

V ^32Re(A) 


-(A-bgi 9 i-f 3292) -V23iIm(A) 0 

—K '/2giKe(A) 0 

0 0 

y/2gilm{A) -w^ill + ci(t)] -71 
0 0 0 

y232lm(A) 0 0 


-V2g2lm{A) 0 \ 

y232Re(A) 0 

0 0 

0 0 

0 UJm2 


-Wm2[l + C2(t)] - 72 / 
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